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Abstract 



Infrared transmission measurements on the 5=1 antiferromagnetic 
chain compound NENP in applied magnetic fields show a sharp absorp- 
tion line at the field-shifted Haldane gap. This violates a wave- vector 
selection rule of the Hamiltonian normally used for NENP, as the gap 
excitations occur at the Brillouin zone boundary. We argue that the 
crystal structure admits terms which can explain the absorption lines. 
In addition, in an applied field, staggered orientations of the g-tensors 
produce a staggered magnetic field. This can explain the observation of 
a finite gap at all applied fields. 
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The one dimensional S — 1 antiferromagnet NENP has provided detailed exper- 
imental evidence |l|]- || for the current theoretical picture [0- [10| of S = 1 anti- 



ferromagnetic chains, and continues to be the subject of a number of investigations. 
However, there remain some puzzling experimental results. One of these involves 
infrared transmission measurements in an applied magnetic field in which a sharp 
absorption line has been observed at the field-shifted Haldane gap. This violates 
the wave- vector selection rule present in the simple model Hamiltonian generally em- 
ployed to describe the system. In addition, the measured gap does not vanish at a 
critical magnetic field but appears to turn around after reaching a finite minimum 
value. We propose that the data can be understood by including in the Hamiltonian 
terms arising from a staggered crystal structure which lift the wavevector selection 
rule. When a uniform field is applied, the staggered g— tensor effectively produces a 
staggered magnetic field, preventing the energy gap from closing. 

We begin by comparing the resonance energies from the transmission experiments 
H with the field-shifted gaps at the zone boundary from neutron scattering exper- 
iments Q (Figure 1(a)). The transmission experiments are at a temperature 1.6iT, 
much smaller than the gap energy (14K), so that the absorption is dominated by 
transitions out of the ground state. The correspondance of the data demonstrates 
that light creates excitations out of the k = ground state into k = n/a excited states 
at the zone boundary. This violates the momentum selection rule of the following 
Hamiltonian, generally employed for NENP: 

H = J £ Si ■ S i+1 + D £(5?) a + E £ [(Sf ) 2 - (S* f] - // B £ ^ • gi? (1) 

i i i i 

where J is the iVi 2+ -iVi 2+ exchange, D and E are single ion anisotropics and g is 
the gyromagnetic tensor. The z— axis is chosen along the chain axis and a is the 
spacing between the iVi 2+ ions. However the correct crystal symmetry of NENP at 



low temperatures is given by the space group Pn2ia ||. This contains the symmetry 
element 2 1; which is a two-fold screw axis consisting of a displacement a along the 
chain, combined with a 180° rotation about the chain. The screw axis arises from the 
alternating orientation of the NO2 groups between the iVz 2+ ions, which is indicated 
schematically in Fig.2. The arrangement of NO2 groups on any given chain leads to 
a uniform electric dipole moment along the chain, and an alternating electric dipole 
perpendicular to the chain. It also implies that the D and the g tensors of two 
successive Ni 2+ atoms on a chain have differently oriented principal axes. The correct 
Hamiltonian should therefore be 

H = Jj2Si-S i+1 + £$-D li S i -// B <S ? i -g 1 #] + [£ ^ ■ D 2 £ - fji B ^ ■ g 2 #] (2) 
» ie(i) ie(2) 

where 1,2 refers to the two sublattices on the chain, corresponding the two orienta- 
tions of the principle axes of the D, g tensors, and we have neglected the staggered 
anisotropy of the exchange interaction. The presence of the twofold screw axis implies 
that D 2 ,g 2 are obtained from Di,gi by a 180° rotation about the chain axis. As 
a function of the applied field, ([]]) and (0) have different phase diagrams: for (fl|), 
there exists a critical field strength H c || for any orientation of the field, such that 
for H < H c there is no spontaneous staggered polarisation. As H crosses H c , there 
is an Ising transition into a state with staggered polarisation. For the Hamiltonian 
the symmetry is explicitly broken, and a staggered polarisation exists for any 
H 7^ 0. Moreover, we shall argue below that consistent with the actual symmetry 
of NENP and Eq.2, there exists a non-zero electric dipole matrix element leading to 
infrared absorption at the Haldane gap even at H = 0, whose magnitude could be 
large enough to explain the observations. Magnetic dipole transitions implied by Eq.2 
could also be important at large H. 

The Zeeman term in (2) can be written as ■ Si + (—l) l h s ■ Si where h = 
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/ieif(gi + g2)/2 is the effective uniform field and h s = psH(gi — g2)/2 is the effective 
staggered magnetic field. h s mixes the ground state with the zone boundary excited 
states and prevents the energy gap from closing at a critical field. This staggered 
field has been invoked previously by Chiba et al. || to understand the field-dependent 
splittings of NMR frequencies observed below the critical field in NENP. These authors 
estimate Sg yz /g = 0.01. 

We have employed an approximate model in order to estimate the quantitative 
effects of the staggered magnetic field on the triplet energies. For simplicity, we ignore 
here the staggering of the D— tensor. We model the low lying excitations via a bosonic 
quantum field theory, which is modified somewhat from the one proposed by Affleck 
H . In absence of a staggered field, the field-shifted gaps from our Lagrangian coincide 
with results of a fermionic field theory proposed by Tsvelik || , which appears to give 
better agreement with experiment than the unmodified boson theory. We do not use 
the fermion theory because there the staggered polarisation is a non-local operator. 
The Lagrangian we use is (h = 1): 

• 2 />; /A, 

€ iji 

+ — e ijk^klm\ j A2 m hihl<pj(j) m ~ A(0 2 ) 2 - h s ■ p(j) , (3) 

ZV ijklm V A fc 

where the spin-wave velocity v and the gaps Aj (i = 1, 2, 3) at the zone boundary 
are taken as phenomenological inputs. We use the values v = 110K, Ai = 15. 7K, 
A 2 = 13. 6K, A 3 = 29K Q, 0. p is a constant relating <fi to the staggered polarisation 
m s = Y,iLi(—]~ySi, so that m s = p<p. We estimate p from the numerically 
determined staggered susceptibility [0 Xs = 20/ J (gps = 1) for the isotropic model. 
In the isotropic limit of the quadratic Lagrangian, \s — p 2 v/A 2 . Using v = 2 J, A = 
0.41 J, we estimate p = 1.3. The above Lagrangian with h s = 0, Ai = A 2 , produces 
the same field dependent gaps at the zone boundary as the fermionic Lagrangian in 
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c= I dt dx \Uf) 2 -1(f) 2 -Y%^T 



v V A, 



0. As an illustration, for h s = 0, h = hie, the lowest magnon branch has energy 
given by 

J(k) = e + h 2 + i(A= + A*) 

At fc = 0, this formula is identical to Eq.l6c in ||, while the formulas differ slightly 
for finite k. We remark that Date and Kindo have discussed a model where the 
magnon energies at the zone boundaries are those of a single 5=1 spin in presence of 
uniaxial anisotropy. After proper identification of the zero-field gaps, the field-shifted 
gaps calculated from such a single- ion Hamiltonian are identical to (§) and Eq.l6c in 

To estimate A, we need to know the non-linear staggered susceptibility. In the 
isotropic limit, the mean field staggered polarisation is m s = Xsh s — A(xsh s ) 3 + O (hg) , 
where Xs — p 2 v/A 2 and A = 4Xv/(A 2 p 2 ). From perturbation theory, we obtain for 
large h s , m s = 1 — (J/h s ) 2 + 0((J/h s ) A ). Matching the large and small h s behaviour 
using a simple functional form m s = Xsh s / 1 {l+Bx 2 s h 2 +xih 4 s ) l ^ A , we obtain an estimate 
A = 1, A = A 2 p 2 /(4v). Assuming A = 0.41J, v = 2J, J = 55K we have a rough 
estimate, A = 3.7K. 

Finally, we need to estimate the staggered field produced by an applied uniform 
field h. Consider the case where h || z. Taking into account the observations of 
Chiba et. ai, we calculate the field dependent gaps for (i) h = hz, h s = 0.01/iy. 
A transverse field applied along the a— axis makes an angle of 32° with the x— axis 
of an individual chain. We therefore calculate field dependent gaps for (ii) h = 
/i(xcos(32°) + ysin(32°)), h s = 0.01 sin(32°)/iz. The mean field values of the zone 
boundary gaps are plotted in Figures l(a),(b), and Figure 3 shows the corresponding 
staggered polarisations. 
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We now consider the interaction Hamiltonian responsible for the absorption. The 
momentum selection rule can in principle be lifted by the presence of localised impu- 
rities or by having chains of finite length. However, in this case, excited states would 
be generated with all possible momenta, producing a threshold instead of a sharp 
line. The observed lines are sufficiently narrow so that this seems to be an unlikely 
explanation. We will restrict our attention to Hamiltonians that preserve the crys- 
tal symmetry. We consider two cases, (a) interaction Hamiltonians that violate the 
translational symmetry of (1), and (b) interaction Hamiltonians that maintain the 
translational symmetry but where the absorption is allowed in presence of a staggered 
polarisation when H ^ 0. In both cases, we consider (i) magnetic and (ii) electric 
dipole transition Hamiltonians. We present below results only for the dc field along 
the chain axis. 

In case (a)(i), the transition Hamiltonian is Hj(t) = —^BH r f{t) ■ 5g ■ l) l S\ 
where H r f is the magnetic field of the infrared radiation. The absorption intensity for 
the two lower branches from this Hamiltonian should be a maximum when H r f || z, 
and should vanish for H r f _L z For H r f \\ z, the absorption coefficient is 

a(u) = ^< S ^b)\ |» 1 J2((-l)^Si(t)Si)e^dt (5) 

where n is the density of nickel ions, (•) is the ground state expectation, and 
Sg yz /g ~ 0.01 ||. At zero applied field, assuming axial symmetry, ((— 1) J_ - J ( 5* (t)S^) ~ 
(l/3)(e~ lA+t / h + e~ %A ~ t ' n )e~ Tt f {\i — j\), where /(0) = 1, and f(i) decays exponen- 
tially for large % with the transverse correlation length £. At resonance, this yields 
for the lowest branch, a ~ (3 x lO~ 5 m~ 1 )Q(!;/a) cos 2 (#m), where Q = A + /hT and 
H r f ■ z = H r f cos(#a/). We estimate the Q value at low fields from the experimental 
data to be be crudely 100, and also £/a « 10 0. This yields a w 3 x lO^m^ 1 . 
Assuming a sample thickness of 1 mm, this gives an absorption dip of 0.003% of 
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transmitted light, small compared to the observed absorption of ~ 0.5%. 

In case (b)(i), the transition Hamiltonian is Hj(t) = —fiBH r f(t) ■ g • J^iS" 1 - The 
orientation dependence of the absorption is the same as in (a)(i). For H r f || z, the 
absorption coefficient can be obtained from Eq.(^) with the replacements 5g yz — > 
g zz , (— l) 1- - 7 — > 1. At low applied fields, noting that J2Si ~ (^/ v )J2 ( Pi x 4>% 
we find a mean field estimate of the absorption at the lowest gap energy a ~ 
«m<5(^/o)(A+/2J) 2 < <p y > 2 cos 2 (#m), where a M = Hon(g zz ii B ) 2 /Tic = 0.3 m~ l With 
parameter values as above, this value of a is again too small to explain the observed 
low field intensities. At high fields, however, the staggered polarisation becomes quite 
large. The lowest mode then predominantly consists of angular fluctuations, and we 
may expect a linearised spin-wave calculation to give qualitatively correct values of a. 
Such a calculation gives a ~ aM<3(2J(</> 3/ )/A + ) cos 2 (9 M ). The line broadens at high 
fields, and taking parameter values to be Q ~ 10, (<fi y ) ~ 0.3, A + w 5K, 6 M = 0, we 
obtain a ~ 200m~ 1 , implying an absorption coefficient of 20% for a sample 1 mm in 
thickness. This is consistent with the experimental absorption at high fields. 

We now consider electric dipole transitions. The iVi 2+ positions lack a center of 
inversion symmetry 0, and electric dipole transitions are therefore allowed within 
the spin manifold. An appropriate Hamiltonian for case (a) (ii) is a time-dependent 



Dzyaloshinskii-Moriya (DM) Hamiltonian [JTT 



Hf ec = P E rf (t) x n perp ■ XX-1)\$ x S i+1 , (6) 

i 

where E r f is the oscillating electric field from the light wave and n perp is a unit vector 
along the staggered component of the internal electric field Ei n . The static DM 
exchange is of the order of J a ~ J8g/g ~ IK. We write p = rjJA/E in , where 
Tj is a dimensionless constant, and we estimate E in as the field produced by one 
electronic charge at a distance of 3A This implies p = rjlO~ 33 C — m. The absorption 



in the lowest branches from the above Hamiltonian vanishes when E r f _L z and is a 
maximum when E r f \\ z. For the latter case, 



where = i x n perf , • (—l) % Si x S i+ i. From the field theoretic description it follows 
that l) l Si x S i+ i ps (l/v)Y,d(f)i/dt. Noting this relation, we estimate a for the 
gap energy to be a ps oieQ(£,I a J )(k. + /2J) 2 cos 2 (#), where E r f ■ z = E r f cos(9e) and 
cue = np 2 /e hc = Arf x 10~ 3 m _1 . Taking typical values at low fields to be Q ps 100, 
A + ps lOfT and £/a ~ 10, we obtain a ~ 0.04r] 2 m _1 . In order to obtain an absorption 
comparable to the experimental value of 0.05% for a 1 mm thick sample, we must 
take rj = 10, which is large, but does not seem impossible. At high fields, a linearised 
spin wave theory yields for the lowest mode a ps a s Q(2J/A + )(0 ?/ ) 2 cos 2 (9 E )- For 
Q ps 10, 7] = 10, ((f) y ) ps 0.3, A + ps 5fT, 6e = 0, we obtain an absorption coefficient of 
10% for a sample 1mm in thickness, comparable to the experimental absorption at 
high fields. 

The Hamiltonian for case (b)(ii) includes a time dependent symmetric exchange 
and a time dependent on-site anisotropy. The absorption from these terms have the 
same selection rules at high fields as the case (a)(ii). 

To summarise our observations on the absorption intensities: while these do not 
strictly vanish for any orientation of the polarisation plane of light, the magnetic 
dipole or electric dipole absorption intensities for a particular line can individually 
vanish for special orientations of the polarisation plane. The magnetic dipole absorp- 
tion at low fields is estimated to be too weak to explain the observations, while the 
spin wave estimate at high field for the lowest branch produces numbers closer to the 
observation. We cannot estimate the electric dipole absorption with accuracy, but our 
crude estimate indicates that it can produce plausible absorption strengths at all the 
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applied fields. A priori, there is no reason why the two absorption mechanisms should 
have the same strength, so that one might expect significant angular dependence in 
the intensities. Careful studies of the angular dependence of the intensities should 
clarify the absorption mechanisms. The strong increase in the absorption intensities 
at high fields is presumably related to the appearence of a large staggered moment. 

In conclusion, the experimentally observed infrared absorption in NENP at low 
temperatures M can be understood after taking into account the crystalline sym- 
metry of the material. The presence of a staggered, anisotropic g-tensor produces a 
staggered internal field when a uniform magnetic field is turned on. This prevents 
the gap from closing at a critical field. Mean field calculations of the gap energies 
qualitatively agree with the experimental observations. Transition hamiltonians can 
be found consistent with the reported crystalline symmetry of NENP that break the 
translational symmetry of the spin chain, and can explain the observed transitions. 

The conclusions of our work are corroborated by recent spin resonance exper- 



iments [E2| where the temperature was varied. The absorptions corresponding to 
those described in || rapidly disappear with increasing temperature, consistent with 
the depopulation of the ground state. As the temperature is raised, absorptions cor- 
responding to flipping the spins of excited magnons increasingly become visible. One 
remaining mystery about the absorption, however, is the splitting of each absorption 
peak observed in the experiments of Ref . || . 
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FIGURES 



In (a), the field shifted gaps from neutron scattering (circles) for h || z are com- 
pared to the positions of the absorption dips in the transmission experiment (crosses) 
. (b) corresponds to field applied transverse to the chain, for which only neutron 
scattering data (circles) is available. The solid lines in (a),(b) result from the mean 
field calculation described in the text. 

Schematic diagram of the backbone of an NENP chain, showing its staggered 
structure. 

Mean field staggered polarisations corresponding to the gaps in Fig.l. The NMR 
lineshifts from Chiba et.al, are shown as dots. The lineshifts have been scaled by a 
constant to make one point lie on the theoretical curve. 
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